Excision boundary conditions for the conformal metric 
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Shibata, Uryti and Friedman recently suggested a new decomposition of Einstein's equations that 
is useful for constructing initial data. In contrast to previous decompositions, the conformal metric 
is no longer treated as a freely-specifiable variable, but rather is determined as a solution to the 
field equations. The new set of freely-specifiable variables includes only time-derivatives of metric 
quantities, which makes this decomposition very attractive for the construction of quasiequilibrium 
solutions. To date, this new formalism has only been used for binary neutron stars. Applications 
involving black holes require new boundary conditions for the conformal metric on the domain 
boundaries. In this paper we demonstrate how these boundary conditions follow naturally from 
the conformal geometry of the boundary surfaces and the inherent gauge freedom of the conformal 
metric. 
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I. INTRODUCTION 

Most numerical relativity applications solve Einstein's 
field equations with the help of a 3-1-1 decomposition, 
which slices the spacetime into a foliation of spatial hy- 
persurfaces and splits the equations into a set of con- 
straint equations and a set of evolution equations P, Q ■ 
The constraint equations restrict the geometry of each 
hypersurface, representing an instant of constant coordi- 
nate time, while the evolution equations determine how 
the geometry changes from one hypersurface to the next. 
In its simplest form, constructing initial data for an evo- 
lution calculation therefore requires finding a solution to 
the constraint equations. 

The Einstein constraints form a set of four equations - 
one in the Hamiltonian constraint and three in the mo- 
mentum (or vector) constraint - and therefore can deter- 
mine only four of the initial-data variables; the remaining 
variables are freely specifiable and have to be chosen inde- 
pendently before the constraint equations can be solved. 
A decomposition of the constraint equations separates 
the freely specifiable variables from the constrained ones. 
Given a particular decomposition, the construction of ini- 
tial data then entails making well-motivated choices for 
the freely-specifiable variables that encode the physical 
characteristics of the system one wishes to model. 

Most decompositions of the constraint equations con- 
formally decompose the spatial metric (see Sec. Ill Al 
for details) and treat the conformally related metric as 
a freely-specifiable variable, meaning that this "back- 
ground" geometry can be chosen arbitrarily. In most 
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applications, the conformally related metric is then cho- 
sen to be fiat. This choice simplifies the equations dra- 
matically, and while it may encode a certain amount of 
physically unrealistic gravitational radiation in the initial 
data - which often manifests itself as "junk radiation" 
when the data are evolved - it is rarely clear a priori how 
to make a better choice. 

An attractive alternative was recently proposed by Shi- 
bata, Uryvi and Friedman ([^; hereafter SUP). In their 
new initial-data decomposition, the conformal metric is 
no longer treated as a freely specifiable variable, and 
is instead determined during the solution of the equa- 
tions. Aside from the trace of the extrinsic curvature, the 
new set of freely specifiable variables includes only time 
derivatives of the metric and extrinsic curvature. For the 
construction of equilibrium or quasiequilibrium data it 
is much more natural to specify the time derivative of 
metric quantities rather than the quantities themselves, 
which makes this new decomposition very appealing. 

So far, this formalism has been used only for bi- 
nary neutron stars When solving the constraint 
equations for black holes, the black hole interior is of- 
ten excised to avoid singularities, which requires suit- 
able boundary conditions on the black hole horizons (see 
[1 i, 0, i, B El [HI)- In the context of constrained 
evolution, excision boundary conditions for all variables, 
including the spatial metric, were explored in Ref. p^ . 
Because the spatial metric is evolved in this case, and 
the spacelike boundaries used have no incoming charac- 
teristics, no boundary conditions on the spatial metric 
were needed. However, since the new SUF formalism 
treats the spatial conformal metric as a constrained vari- 
able rather than as either a freely-specifiable metric or an 
evolved quantity - and hence as a solution of an elliptic 
equation - we need to provide suitable boundary con- 
ditions for the conformal metric on black-hole excision 
boundaries before this formalism can be applied to black 
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holes. In this paper we demonstrate how these boundary 
conditions can be formulated quite naturally in terms of 
the conformal geometry of the excision boundaries and 
the inherent gauge freedom in the conformally related 
metric. 

We found it useful to develop this formalism in terms 
of a reference metric approach. While this approach 
is not new to general relativity, it has not been used 
widely within the numerical relativity community. We 
will therefore present a full description of the reference 
metric approach and write all of the Einstein evolution 
and constraint equations in terms of this formalism. 

Finally, we must keep in mind that for coupled non- 
linear elliptic equations, it is not always clear whether 
or not a given set of boundary conditions are indepen- 
dent and lead to a well-posed elliptic system. In Refs. Q 
and 0] , two different physically motivated boundary con- 
ditions were developed for the lapse for use with the ex- 
tended conformal thin-sandwich equations (see Sec. lIIB| ). 
When combined with the boundary conditions for the 
remaining elliptic variables, however, neither condition 
gave rise to an independent set of boundary conditions, 
and hence left the system degenerate (see Refs. 6] and 
[m , and also Ref . [l^ for a general discussion of elliptic 
systems). While we will derive a well-motivated set of 
boundary conditions for the conformal metric, we have 
not yet shown that they lead to a well-posed elliptic sys- 
tem. 

This article is organized as follows. We provide an 
overview of the problem in Sec. |TT] and derive the ref- 
erence metric approach in Sec. IIIII We then derive the 
boundary conditions on the conformally related metric 
in Sec. IIVI and briefly summarize in Sec. |Vl We also 
include a complete list of Einstein's equations in the ref- 
erence metric approach in App. |^ 



II. OVERVIEW 

A. 3-1-1 and conformal decompositions 

We begin with the spacetime metric written in the gen- 
eral 3 + 1 form 

ds2 = -a'^dt^ + 7y (ds* + f3'dt){dx^ + (3^dt), (1) 

which treats the full spacetime as a foliation of spacelike 
slices with timelike unit normal vector n'^. In Eq. ([1]), 
a is the scalar lapse of time, /3' is the shift vector, and 
7ij is the metric of a spatial slice. In terms of this de- 
composition, the vacuum Einstein equations are written 
as 

dtK,, = -V,Vja + a[R,,-2KuK^j+KK,j] (2) 
^a^J = -2aK,j + 2%l3j^ (3) 



= R + - K,jK'^ 
= Vj{K'^ -f^K) 



(4) 
(5) 



Equations Q and ([5]) are, respectively, the Hamiltonian 
and momentum constraint equations, and Eq. Q is the 
evolution equation for the extrinsic curvature Kij which 
is itself defined by Eq. ^ . The covariant derivative com- 
patible with the spatial metric is written as V^, the 
Ricci tensor and scalar curvature as Rij and R, and dt is 
the derivative along the time vector 



(6) 



Finally, it is often convenient to decompose the extrinsic 
curvature into its trace K = ^"^^ K^j and trace-free parts 



(7) 



It is also useful to conformally decompose the spatial 
metric and other quantities. In particular, we will make 
use of the conformal factor ip which allows us to define 
the conformal metric 7ij, conformal trace-free extrinsic 
curvature Aij and conformal lapse a via 



7ij — ^ ^ij 

Aij = ip-'^A,_ 



(8) 
(9) 
(10) 



In terms of these conformal variables and the covariant 
derivative Vi, Ricci tensor Rij and scalar curvature R 
compatible with the conformal metric jij, we can write 
the various constraint and evolution equations as 






dtA'' 



WjA'^ - li^^f^VjK 
~2aAj + {t[4,-^|3])^J 
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%l;^a{R}^ - \f^R) 



M,A^' (11) 
(12) 
(13) 
(14) 



(0«5) 



(V(* lni/;)V^') ln(V'^(5) 

_i^»J(VfelnV;)Vfeln(i/'^5) 



dtK = -V'~^V2(V'^a) -f |V^a^-K ^V^aiiT^ (15) 



dti^ = i7/;(Vfe/3'=-|-6/3'=Vfclni/;-7/;65if). 



(16) 



Here the conformal longitudinal derivative of a 1-form 
(L \j is defined by 



(17) 



and we note that the term -i/;"^ appears in its argument 
in Eq. (fT5|) because the index of the shift /3* has to be 
lowered with the physical spatial metric 7y rather than 
with the conformal metric 7y . 
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B. Current initial data formulations 



The most commonly used decompositions of the con- 
straint equations are the conformal transverse-traceless 
and the conformal thin-sandwich decompositions (see, 
e-g-, [3, d, EBl for reviews). 

The conformal transverse-traceless decomposition 
solves the Hamiltonian and momentum constraints, 
Eqs. ([TT|) and ([H]), and fixes the conformal metric 7^- as 
well as the mean curvature K and the transverse-traceless 
parts of Aij as the freely-specifiable variables. If the con- 
formal metric is chosen to be flat, 7^ = /y, and the 
initial hypersurface to be maximal, so that K = Q, the 
momentum constraints decouple from the Hamiltonian 
constraint and become linear. Vacuum solutions describ- 
ing one or multiple black holes, known as Bowen-York 
solutions [Tt'I, can then be found analytically and form 
the basis for both _2uncture [H, [l^ Ho] and conformal- 
imaging solutions HH to the Hamiltonian constraint. 

In the conformal thin-sandwich decomposition [23| the 
evolution equation for the conformal metric, Eq. is 
used to replace the transverse-traceless parts of Aij with 
the time derivative of the conformally related metric in 
the set of freely-specifiable variables. Effectively, we can 
solve Eq. (|13p for Aij and then insert this equation into 
every occurrence of Aij in the other equations. We now 
fix the trace of the extrinsic curvature K, the conformal 
metric 7ij, its time derivative dtjij, as well as the con- 
formal lapse a as the freely specifiable data and solve 
the Hamiltonian constraint (fTT|) for the conformal factor 
ip and the momentum constraint (fT^ . via (|13p . for the 
shift p\ 

The extended conformal thin-sandwich decomposition 
[2^ incorporates the time derivative of the trace of the 
extrinsic curvature, Eq. (fT5l) . If we take dtK as specified, 
this equation yields an elliptic equation for the conformal 
lapse. Essentially, we now replace the conformal lapse 
with dtK in the subset of the freely specifiable variables, 
and we obtain one additional elliptic equation which must 
now be solved for the conformal lapse. While this change 
of perspective results in an additional elliptic equation 
which must be solved in order to construct initial data, 
this added burden is offset by the fact that the set of 
freely specifiable data now takes on a particularly at- 
tractive form. Specifically, we now specify the conformal 
metric and its time derivative, along with the trace of 
the extrinsic curvature and its time derivative. In many 
situations, in particular as mentioned before for the con- 
struction of equilibrium or quasiequilibrium initial data, 
it is advantageous to specify the time derivative of some 
quantity rather than some other field. (See [23 . [25l .[26l.[27j 
for a discussion of the uniqueness issue in the extended 
conformal thin-sandwich formalism.) 



C. Determining the conformal metric 

In all of the initial-data methods sketched out in 
Sec. IIIBl the conformal metric is taken as part of the 
subset of freely-specifiable initial data. Roughly speak- 
ing, the five degrees of freedom that are fixed by spec- 
ifying the conformal metric include the initial choice of 
the spatial gauge and two dynamical degrees of freedom. 
Therefore, fixing the conformal metric strongly affects 
the initial gravitational radiation content of the initial 
data. 

In most cases, the conformal metric is chosen to be 
flat. But even when it is not chosen to be flat, in all but 
the most trivial cases, the chosen metric is not completely 
compatible with the physics that one wishes to build into 
the initial data. The result is that undesired "junk" grav- 
itational radiation is built into the initial data. While this 
defect in the data is usually small and has little effect 
on the gross physics one wishes to simulate, it can have 
a significant effect on detailed comparisons with post- 
Newtonian methods [1^ , and may impact future param- 
eter estimation efforts. 

It is desirable to find a way to construct initial data in 
which the conformal metric is not fixed a priori, but is 
constructed in a way that is consistent with the physics 
one wishes to simulate and that eliminates, or at least 
reduces, the junk radiation. Following in the spirit of 
the extended conformal-thin-sandwich approach, SUF 
noticed that the equation for the time derivative of the 
trace-free part of the extrinsic curvature, Eq. (jl4p . can be 
written as an elliptic equation for the conformal metric if 
the spatial gauge is imposed in a suitable way and if we 
take dtA'^^ as freely-specified data. In close analogy to 
going from the conformal transverse-traceless to the ex- 
tended conformal thin-sandwich decomposition, namely 
using the evolution equation for the conformal metric to 
replace Aij with dt"fij, and the evolution equation for 
the trace of the extrinsic curvature to replace a with 
dtK as freely specifiable variables, we now use the evo- 
lution equation for the extrinsic curvature to replace 7^ 
with dtAij. This formalism is very attractive since it is 
again more natural to make choices for time derivatives 
of functions than for the functions themselves. 

To date, this approach for constructing initial data has 
been implemented successfully only for binary neutron 
stars ^] . The main goal of this paper is to develop the 
formalism necessary for applications to black-hole initial 
data where the black-hole interiors are excised from the 
computational domain. When black-hole initial data are 
computed using excision methods, boundary conditions 
for all of the initial data that are determined by elliptic 
equations must be applied on the excision boundaries. 
For the case of the extended conformal-thin-sandwich 
initial data, the boundary conditions have been worked 
out and thoroughly tested 5, 6, 7, 8]. In the context of 
the new formalism of SUF, however, additional boundary 
conditions are required for the conformal metric. 
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III. REFERENCE METRIC APPROACH 

Our goal is to derive a system of equations that will de- 
termine the conformal metric together with some other 
quantities. Many of the operators that need to be in- 
verted to construct the conformal metric depend on the 
conformal metric themselves. It is therefore convenient 
(though certainly not necessary) to use a reference met- 
ric approach wherein, in addition to the conformal metric 
7ij , we also associate some appropriate fixed metric with 
our solution manifold. We can then formulate the oper- 
ators in terms of this fixed metric, which greatly simpli- 
fies the inversion of the operators. This approach is not 
new (cf Refs. [H, HODj but is not widely used within the 
numerical relativity community. A notable exception is 
[31], whose formalism shares many elements with ours. 
We will give a basic outline of the approach in this Sec- 
tion, and hst the complete set of Einstein's equations in 
Appendix [X] 



The latter can be rewritten as 
SRij — —iti-j — g^™'gk{iRj)em'' — -^9^"^^e^mgij (23) 

+ V(, {gj)kg''"'sr\^^) ~ gpkg'"'sr^raSrp,, 

+ ^/™5"''{(Vj5mp)V£.gj„ -I- (Vj.9mp)Vf5„ 
+ (y tgin)^ mgjp — tgm)^ pgjm 

As it turns out, this form is particularly useful because 
all second derivatives of the gij have now been absorbed 
into only two terms. The first of these two terms, 
f/^™ V^V^ (7y , forms an elliptic operator acting on gij as 
long as both gij and gij are sufficiently well behaved, 
while the second of the two terms, V(j {gj)kg^"^Sr''em) , 
can be eliminated by virtue of a suitable gauge choice, as 
we will discuss in the next section. 



B. Reference metric and gauge choice 



Basic outline 



We assume that our initial data hypersurface is rep- 
resented at a basic level by a manifold with coordinates 
and coordinate maps defined everywhere. We associate 
with this manifold two metrics which are not necessarily 
the same: gij and gij. Each metric has an inverse and 
covariant derivative such that 



gjkg^'^S} and g.kg"' ^ (18) 
Vfcgy = and Vkfjij = 0. (19) 



The difference between two connections is a tensor which 
can be written as 



1 



(20) 



-,9 



ki 



Vi.9j£ + Vjgie - Vegij 



The differences between the Riemann and Ricci tensors 
for each metric can be written as 



SR 



■ijk 



^ij k ^ijk 



(21) 



■ ST^kiST^mj — ST"^kjST^mi 



Riij — Riij ' 
I 



(22) 



We will assume that we are dealing with one or more 
black holes in an asymptotically flat initial-data hyper- 
surface. In this case, the solution domain is topologically 
with a "ball" cut out for each excised black hole inte- 
rior. It is therefore appropriate to take the reference met- 
ric (jij to be a flat metric (not necessarily in Cartesian 
coordinate) . For consistency with our previous notation, 
the "other" metric gij will be the conformal metric we 
wish to determine, 7ij . To make the notation as clear 
as possible, we will denote the covariant derivative com- 
patible with the flat reference metric as and so the 
difference of connections becomes 



^ ij 

[V^loi^+VJl^l~Vll^o] 



(24) 



We next impose a spatial gauge condition by setting 
the contractions of the connection coefficients equal to 
some predetermined gauge source functions ^ 



5f^ = f^5f^ij =-- 



1 



\/det7 



=V,{V^l''')^V\ (25) 



For convenience and ease of reading, we will denote 
the determinant of the conformal metric as either det 7 
or 7. One possible choice is = 0, in which case 
we would obtain a "spatial conformal harmonic gauge" , 
but we will leave arbitrary for generality (com- 
pare also the "generalized Dirac gauge" of fsill). With 
the gauge source functions V'^ given as specified func- 
tions of the coordinates, we see that the only remaining 
second-order term in Eq. ((23l) is the first one, bring- 
ing the conformal Ricci tensor into a familiar elliptic 
form. This or similar properties of the Ricci tensor have 
often been utilized before, both in a four-dimensional 
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(e.g. [3^. [bsI . [3I. 35 J36 . STLISSj l and a three-dimensional 
context (e.g. 30. l3ll. l39l. l40j|). An attractive feature of 
the reference metric approach is the fact that differences 
of connection coefficients ST'^ij are tensors, so that the 
gauge source functions V'^ become vectors and take a 
gauge-invariant meaning. 

In our case, the conformal Ricci tensor and conformal 
Ricci scalar take the form 



(26) 
(27) 



1 ~, 
2 



%e [C 



(compare Eqs. (44) - (47) in [3l|). Here fA is a general- 
ized Laplacian defined by 



(28) 



and the symmetric tensors B^^ , C'-' , I?*-' , and are 
quadratic combinations of first-derivatives of the confor- 
mal metric defined by 



B'^ EE 7mn7"(yfc7'"')y£7^" = 

EE 7m„(yfc7™<07^'^V£7'" = 



(29) 
(30) 
V'' (31) 
(32) 



We note that the derivatives of the gauge source func- 
tions have been written in terms of the conformal covari- 
ant derivative Vfc, and not the flat covariant derivative 
^fe. Also, in the future, we plan to solve directly for 
the inverse conformal metric 7'-^ and so the conformal 
Ricci tensor and the various derivatives of the metric are 
written with the indices raised (compare [3l|). We list 
the complete set of Einstein's equations in the reference 
metric form in Appendix \K\ 



IV. BOUNDARY CONDITIONS 

Our goal in this section is to determine the proper 
boundary conditions for the inverse conformal metric 7*-' 
required on the excision boundaries when we solve an 
elliptic equation for this variable. We begin with some 
simple counting arguments. 

The conformal metric 7ij is a member of a conformal 
equivalence class of metrics which each have five indepen- 
dent degrees of freedom. To fix upon a particular member 
of this class, we can fix the determinant of the conformal 
metric det(7). As mentioned previously, these five de- 
grees of freedom are usually thought of as containing the 
two dynamical degrees of freedom and the three spatial 
gauge degrees of freedom in the metric. Of course, this 
splitting is not clean and, in general, they mix with the 
sixth constrained (or longitudinal) degree of freedom of 
the full spatial metric 7^^. However, this rough splitting 
suggests how we should fix the boundary conditions. 



First, the excision surface is a closed 2-surface with 
topology. The metric on any such surface is conformally 
equivalent to the unit sphere. Since the metric induced 
on the excision surface by the physical metric 7^ will 
be conformally equivalent to a spherical metric, it seems 
that we should be able to demand that the metric induced 
on the excision surface by the conformal metric 7^ be a 
spherical metric. In fact, this can be generalized to the 
statement that we have the freedom to fix the metric 
induced on the excision surface by the conformal metric 
7y to that of any metric with S*^ topology. A 2-metric has 
three degrees of freedom, but because of the conformal 
equivalence of all 2-metrics, specifying the induced metric 
really only fixes two degrees of freedom. 

Having fixed two of the five degrees of freedom of 7^ 
on the excision surface leaves three more boundary condi- 
tions that must be specified. Since the conformal metric 
contains three gauge degrees of freedom, it seems natural 
that these three remaining boundary conditions should 
result from the gauge conditions 



A. 2 + 1 decomposition 

To rigorously define and clearly understand the new 
boundary conditions, it is most convenient to rewrite 
the conformal metric in terms of a 2 -|- 1 decomposition 
adapted to the excision surface. We define on the solution 
manifold a foliation of concentric surfaces with topology 
S"^ in the neighborhood of the excision boundary. Each 
surface is defined as a level surface of a scalar function 
r{x^) such that r — tq defines the excision surface. Inde- 
pendently of which metric we are using, we have a normal 
1-form for the excision surface given by 



dir{x^ ) I 



(33) 



In terms of each metric (7^ and fij) we can normalize 
this 1-form and define a unit-normal vector and unit- 
normal 1-form. In terms of the unit-normals to the exci- 
sion surface, we can then define projection operators as- 
sociated with each metric and a corresponding induced 
metric on the surface defined by r — rg. For the con- 
formal metric 7ij, we will denote its induced metric on 
the excision surface as Hab, and for the flat metric /y 
we will denote the induced metric as Sab- Here and in 
the following, upper-case Latin indices denote adapted 
coordinates in the excision surface. 



1. Conformal metric 

Focusing first on the 2-1-1 decomposition of the confor- 
mal metric we write the three-dimensional line interval 



d/ 



A^dr^ -t- liAB (dx^ + A^dr) [Ax^ + A^dr) , (34) 
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and note that the determinant is given by 

det 7 = det ft.. (35) 

Evidently, A and A"^ play the role of the lapse a and shift 
in the more familiar 3+1 decompositions. In "matrix" 
form we write the conformal metric as 



7^j 



hAB 



(36) 



where Xa — hAsX^ , and its inverse as 
A-2 -A-^A^ 
-A-2A^ ft^^ + A-2A^A^ 



(37) 

(38) 
(39) 
(40) 



The unit-normal 1-form is given by 

5, = A a,r(x^) =A[1,0], 

and the unit-normal vector by 

=fKsj = A-i[l,-A^]. 

Similarly, the "radial" vector is given by 

r' = AS' + A' = [1,0], 

where the 3- vector A* = [0, A"^]. 

The induced metric on the excision surface is given by 

hi;j = - SjSj, (41) 

and we define the extrinsic curvature of the excision sur- 
face as 

H,, = hfh'^%Sey (42) 

The covariant derivative compatible with the induced 
metric Hab will be denoted Dq- Because we are us- 
ing adapted coordinates, the Lie derivative in the radial 
direction Cp is equivalent to the partial derivative in the 
radial direction, 



d_ 

dr 



(43) 



Finally, we note that 

drhAB = 2KHab + Da>^b+ DbXa, 



(44) 

which follows from the definition of the extrinsic curva- 
ture. 

Since our gauge conditions (j25p are defined in terms 
of the difference between the conformal and flat connec- 
tions, it is useful to list the connections associated with 
the conformal metric in terms of the 2 + 1 variables. 



!_ rr 
^^rA 
^^AB 



i rB 



9^ In A + A-^Da In A - ^A^A^-ffAi3,(45a) 



1 



DaIuA- ^X^Hab, 



A 



-jHab, 



drX^ - A^t)^ In A + 2AX^M 



X^ D bX"^ — A'^f '^rr; 



pA 



BC 



-A^r 

.A-4f 



BC 



Ji^^DbXc, 



(45b) 

(45c) 
(45d) 

(45e) 
(45f) 



where T^bc is the usual connection constructed from the 
induced conformal metric Hab- 



2. Flat metric 

The 2 + 1 decomposition of the flat metric fij follows 
precisely the same form as that for the conformal metric 
as outlined above in Sec. lIV ATI We simply need to define 
new variables for the analogs of A, A"*, and Hab- While 
it is possible to proceed in full generality, we will restrict 
ourselves to the case where we take the flat metric analog 
of A ^ 1 and A"* ^ so that we have 



and 



1 
Sab 



1 

S^^ 



(46) 



(47) 



Specifying the induced metric Sab is sufficient to fix the 
coordinate transformations between the coordinates on 
the 3-dimensional initial-data manifold and the 2 + 1 
adapted coordinates. For example, the simplest choice is 
to let Sab be the standard metric for a sphere of radius 
r. Then, if the initial-data manifold uses standard Carte- 
sian coordinates, the coordinate transformations are well 
known. 

To proceed, we will need to define the covariant deriva- 
tive compatible with Sab ■ We will denote this derivative 
as ^A- Finally, we can write the connections associated 
with our flat metric as 



pr 


= 0, 


(48a) 


V 

fL rA 


= 0, 


(48b) 


F^AB 


= —-^drSAB, 


(48c) 


pA 

fj- rf 


0, 


(48d) 


pA 

A rB 


- is^^'drSBC, 


(48e) 


pA 

BC 


— si BC: 


(48f) 



where J'^bc is the usual connection constructed from 
the induced flat metric Sab- Note that Eqs. (|48c|l and 
(|48e)) follow from Eqs. (|45c|) and (|45e)) by using Eq. (|44)) . 



B. Gauge boundary conditions 

All the pieces are now in place to fully understand and 
rigorously define excision boundary conditions that can 
be used for the conformal metric. Recall from our discus- 
sion at the beginning of Sec. lIVI that to fix upon a unique 
member from the conformal equivalence class of confor- 
mal metrics jij, we may specify its determinant det(7) 
everywhere in the initial-data manifold. We choose to 
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demand that the determinants of the conformal metric 
and the reference metric are equal. In this paper, we are 
restricting ourselves to the case that the reference metric 
is flat, so we demand 



det 7 = det /. 



(49) 



We also point the reader to Sec. II. B of Ref. [31[ for a 
more in-depth discussion of this point. 

Also recall that Hab and Sab, each having topology 
S"^, are both related to the standard metric of a unit- 
sphere by a conformal transformation and an appropriate 
coordinate transformation. Now, since we are expressing 
both hAB and Sab in the same coordinates, we can de- 
mand that they are related to each other simply by a 
conformal transformation: 



hAB = n\x')SAB. 



(50) 



It follows immediately from Eq. ((35|) . its analog for the 
induced flat metric that det(/) = det(S'), and Eq. (IST 
that 



(51) 



Notice that because of Eqs. (|50l) and the 2 + 1 

version of the conformal metric can be expressed in terms 
of A, A-^, and Sab'- 



ds^ = A^dr^ 



^SABidx^ + \^Ar){dx^ +\^dr). (52) 



These variables hold six degrees of freedom. But recall 
that they incorporate the restriction on the conformal 
equivalence class given in Eq. , so that there are really 
only five degrees of freedom being fixed by the choice of 
these variables. 

Because all 2-metrics are conformally equivalent, spec- 
ifying Sab can now clearly be seen as simply specifying 
the coordinates on the excision surface. The variables 
A and A'^ then determine how these coordinates propa- 
gate off of this boundary surface into the full initial-data 
hypersurface. These three degrees of freedom must be 
chosen so that the coordinates satisfy the gauge condi- 
tions in Eq. (|25p. In terms of our adapted coordinates, 
and making use of Eqs. (I45p and these three gauge 
conditions are 



= --^a.(v/7A-2)+y^(A-2A^) 
+ KdrlnyfS +]-K-^^\'^\'^drSAB, 



(53) 



V7 

A-'\^S^^drSBC, 



(54) 



and again for compactness we define det(/) = / and 
det(5') = S. 

As an explicit example, assume we are using Cartesian 
coordinates on the spatial hypersurface. For simplicity. 



consider the case where the gauge source functions are 
chosen to vanish {V^ — 0) and the excision boundary is 
chosen to be a coordinate sphere of radius r. Then the 
coordinates adapted to the excision surface are standard 
spherical coordinates, we have 



Sab 



r2 
sin^ i 



(55) 



and det(/) — det(S') r sin 6, and the coordinate 
transformations between the three-dimensional Cartesian 
coordinates and our adapted coordinates are well known. 
The three gauge conditions then take the form 



'a,A = 



-{i-A')-UyAX^ + -XAX^] (56) 
r 2 \ r ' 



A^^A - ^ ( 1 + A^ + ^AsA^ ) A-^ (57) 



2 

rA 



A 



A^A^T^ 



BC- 



Note that these conditions are identically satisfied if the 
conformal metric is flat so that A = 1 and A"^ = 0. 

In this example, we impose boundary conditions on the 
six Cartesian components of the conformal metric 7^ as 
follows. First we use standard Cartesian- to-spherical co- 
ordinate transformations to construct the three compo- 
nents of JAB = hAB- In a similar fashion, we construct 
the two components of A^ and the scalar A. The vari- 
ables fiAB, X^, and A are then explicit functions of the six 
Cartesian components of the unknown conformal metric 
7y . We then demand that hAB = Sab /A which con- 
stitute three Dirichlet boundary conditions. Finally, we 
impose the three normal-derivative boundary conditions 
given in Eqs. ([55]) and (|57p . Together, these form a cou- 
pled set of six equations for the six components of 7^ . 
However, we emphasize that the choice of Sab guaran- 
tees that the global condition of det (7) is satisfied, and so 
we are really only constraining the five remaining degrees 
of freedom in the conformal metric. 

Finally, we note that when implementing boundary 
conditions on the Cartesian components of the confor- 
mal metric, one would not actually use the boundary 
conditions as written in Eqs. ([55)1 and (|57p . Instead, one 
would use the gauge conditions of Eq. (P5)l directly in 
terms of Cartesian components. The 2-1-1 decomposition 
we have derived is actually most useful in implementing 
the Dirichlet conditions, and in better understanding the 
nature of the boundary conditions. 



V. SUMMARY AND DISCUSSION 

In this paper we derive boundary conditions for the 
conformal metric that can be applied on black-hole exci- 
sion boundaries. These boundary conditions are needed 
in the context of a new initial-data decomposition re- 
cently proposed by SUF, which provides an equation for 
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the conformal metric rather than treating it as a freely 
specifiable variable. This seems very attractive, since it 
avoids the need to make ad-hoc choices for the conformal 
metric ~ like conformal flatness - and instead computes 
the conformal metric as part of the solution. 

Whenever black-hole initial data are constructed us- 
ing excision methods, boundary conditions for all con- 
strained data are required on the excision boundaries. 
In the conformal thin-sandwich decomposition, bound- 
ary conditions for the conformal factor tp and shift /3* 
had a strong physical motivation that can be naturally 
viewed within the isolated-horizons framework @, 
Interestingly, this is not the case for the boundary con- 
ditions for the conformal metric. 

It might seem that different choices for the induced 
metric Sab should lead to physically different initial 
data. While we have emphasized the conformal equiva- 
lence of all closed 2-metrics, we must recognize that this is 
a statement about the intrinsic geometry of the 2-surface 
only. We must also consider the effect of the choice of the 
extrinsic curvature. Although we are restricted to choos- 
ing a boundary 2-surface with topology 5'^, we have a 
wide range of choices for how this surface is embedded 
in the reference space. For example, we could choose the 
boundary to be either a coordinate sphere or a coordi- 
nate ellipse in a flat Cartesian space, each of which has a 
different extrinsic curvature. Note also that we can make 
this choice independently of the choice of the metric on 
the excision surface; for example, we could choose Sab 
to be the metric of a unit sphere even on a coordinate 
ellipse. Both the choice for the shape of the 2-surface and 
for Sab will affect the extrinsic curvature of the 2-surface 
as defined by Eq. But, it is important to remember 

that these fix the conformal extrinsic curvature of the 
2-surface, and not its physical extrinsic curvature. 

The 2-surface's influence on the dynamical degrees of 
freedom of our initial data are embodied in its physical 
extrinsic curvature via the sheer and expansion of the 
family of outgoing null rays passing through the bound- 
ary 2-surface. The physically motivated boundary con- 
ditions on the conformal factor il) and shift mentioned 
above are, in fact, obtained by demanding the sheer and 
expansion vanish on the boundary. These choices will 
certainly have an affect on the physical extrinsic curva- 
ture of the excision surface. Furthermore, the choice of 
the shape of the excision surface may have some affect 
on the dynamical degrees of freedom of the initial data, 
although the choice for Sab should not. 

As we demonstrate in this paper, it is natural to write 



the conformal metric on the excision surface in a 2-1-1 
decomposition. The induced conformal metric on the ex- 
cision surface must be related by suitable coordinate and 
conformal transformations to the metric on a unit sphere 
- making this choice therefore does not impose any phys- 
ical restriction on the conformal metric. The remaining 
degrees of freedom of the conformal metric on the ex- 
cision surface can be expressed as a radial "lapse" and 
"shift" , and can be determined from the gauge conditions 
that are imposed on the conformal metric. We therefore 
conclude that these boundary conditions affect only the 
gauge degrees of freedom of the conformal metric and in 
no way restrict the dynamical degrees of freedom of our 
initial-data solution. 

While the investigation of these boundary conditions 
was motivated by the need to provide excision bound- 
ary conditions for black-hole initial data, the fact that 
only general geometrical and gauge considerations are 
involved means that these boundary conditions can be 
applied to any topologically spherical boundary. These 
boundary conditions need not be applied at a black-hole 
horizon, and in particular these conditions can be applied 
without approximation on an outer boundary. 

Finally, we repeat our warning from Sec. [J that the set 
of boundary conditions we have derived have not yet been 
shown to be independent and give rise to a well-posed el- 
liptic system. Because the set of elliptic equations are 
non-linear and coupled, it is very difficult (if even pos- 
sible) to determine this analytically. Future numerical 
implementations will clarify this issue. 
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APPENDIX A: EINSTEIN'S EQUATIONS IN REFERENCE METRIC FORM 



For completeness, we include below the conformally decomposed 3-1-1 version of Einstein's equations, expressed in 
a reference metric form, along with a few additional useful equations. In writing these equations, we assume that the 



fundamental independent variables are 

V', <5, (3\ and f\ (Al) 
while we take the following quantities to be freely specifiable 

dtf^, dtA'^, dtK, K, det7, and V\ (A2) 

The conformal metric 'jij and the trace-free conformal extrinsic curvature A^^ are auxiliary variables defined to 
simplify the equations. The conformal metric is derived from the conformal inverse metric 7*^ in the usual way, and 
the trace- free conformal extrinsic curvature is defined by Eq. (|A5|) below. Note that, as in the main text, derivatives 
of the gauge source functions are taken using the conformal covariant derivative whereas derivatives of the 
main variables are taken using the flat covariant derivative y,. We emphasize again that we fix the reference metric 
to be flat (although this can be generalized if needed) and therefore the Ricci tensor and scalar associated with the 
reference metric vanish. 

The Hamiltonian constraint in Eq. (jlip is an elliptic equation for the conformal factor and takes the form 

- ^i^%j{C' - IB'') - j^V'^' + l^-^^jikiA^'^A^' = FVzV' + l^^^V' - 3^VfAln7, (A3) 
and the momentum constraint in Eq. (jl2p becomes 

y^.iy _ ^ ^_^krn^^^~^^~^3JJ^~^n _ A^'^^^Vif' = In 7- (A4) 

The trace-free conformal extrinsic curvature is defined in terms of the time derivative of the inverse-conformal 
metric 

dtf' = 2aA'' - 27'=(Vfe/3^'^ + ll'^Wkl^^ + P^Vkl'^ + \l"' I^^Vu ln7 = -u'^ ■ (A5) 

Note the definition of made for convenience below. The evolution equation for the trace-free conformal extrinsic 
curvature Eq. (jl4p becomes an elliptic equation for the inverse conformal metric: 

dtA^' = [Af' - B'^ - C' - If'jkeC"' + 2£'' - \{V' - lf^%,v''')] (A6) 

- (f 7^' - ^f'j''')Wkmi^'a) ~ (7^(V^7^'^™ - if ™y,7''^)y™(V''5) 
+ 8^«a(7'=(*f - if^7")(yfelnV')y£ln(?/'^a) 

+ + A^'^^kP" + limA'^^'l'^'^VkP' - i"(Vfc/3^') - /3'=7^™i'(Vfc7^>' - 7feM'**"''^' 

while the evolution equation for the trace of the extrinsic curvature Eq. (|15p is an elliptic equation for the conformal 
lapse 

dtK = -^-5^.A(V/5) + ^i>^a%{C' - ^B'') + ^V'^if^ + ZV"'5%-7fe,i'"i^'^ + P'W^K (A7) 
+ V^VVi(V'^a) + lij'^cNiV' - j^i/;25fAln7. 

Finally, the evolution equation for the conformal factor becomes 

dtijj = + 6/3V» In - V^<5X) + ^ij!3% In 7. (A8) 

If we replace A'^ in the momentum constraint (|A4p by its definition in Eq. (|A5[) . we obtain an elliptic equation for 
the shift: 

AP' + ll'^VjVkP'' + i(25i"'")7M7mnf ^^^,7'" - {2aA''')^,uVif' - {2aA''W, Ina - |a^6~^,^^.^ _ y^. ^^9) 
= ~ IVy^P' - P'V.V' - \V'P'Vj ln7 - WaA^'W, ln7 + \l''{VjP''Wk ln7 " \f'{WkP''Wj ln7 

Hl^'iVkP^Wj ln7 - hP'^i^kl^'W, ln7 - If'P^'VjVk ln7 " \n^Wk ln7- 

In all of these equations terms of the form y,; In 7 appear. These terms are defined implicitly by 
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Note that both det(7) and det(/) transform as scalar densities of weight 2, so their ratio is a simple scalar. If we 
choose det(7) = det(/) as we have in defining the boundary boundary conditions in Sec. IIVI then we find that all 
terms involving derivatives of det(7) vanish. 

Finally, we note that the Arnowit-Deser-Misner energy takes on the familiar form: 

Eadm ^-1^ (f Wri^d^S' + T^ (f {V'- V^ In 7) d^ 5, . (All) 
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